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Superconductors with low superfluid density can be described by XY models. In such models 
the scale of the transition temperature T c is largely set by the zero temperature phase stiffness 
(helicity modulus), a long- wavelength property of the system: T c = AT(0). However, the constant 
A is a non-universal number, depending on dimensionality and the degree of inhomogeneity. In this 
Letter, we discuss strategies for maximizing A for 2D XY models, that is, how to maximize the 
transition temperature with respect to the zero temperature, long wavelength properties. We find 
that a framework type of inhomogeneity can increase the transition temperature significantly. For 
comparison, we present similar results for Ising models. 



Many strongly correlated models exhibit either local 
inhomogeneity (whether ordered or disordered) or out- 
right phase separation, and there is experimental evi- 
dence that some degree of local electronic inhomogene- 
ity takes place in various parts of the phase diagrams 
of transition metal oxides such as nickelates, cuprates, 
and manganites.[17| It is important to understand how 
the macroscopic properties emerge out of the mesoscale 
structure, and whether it has detectable consequences for 
the observed phases, such as the technologically impor- 
tant superconductivity observed in some strongly corre- 
lated systems. That is, does local inhomogeneity help or 
harm superconductivity, or is it a side issue entirely? 

For superconductors with low superfluid density, the 
transition temperature is dominated by phase fluctua- 
tions of the superconducting order parameter, and the 
transition may be captured by an XY model. Although 
the transition temperature in XY models is largely set by 
zero temperature, long wavelength properties of the sys- 
tem, dimensionality and inhomogeneity also play a role. 
That is, T c = AT(T = 0) where T is the phase stiffness or 
helicity modulus (proportional to the superfluid density 
in a superconductor), and A is a non-universal number 
of order 1. We focus here on how inhomogeneity may be 
used to maximize A, and therefore enhance the transition 
temperature with respect to the zero temperature, long 
wavelength properties of the system, as compared to the 
uniform case. 

It has generally been expected that inhomogeneity 
should decrease T c , especially to the extent that it in- 
troduces disorder or competing orders. However, this 
intuition has been violated even in conventional super- 
conductors when mesoscale structures were introduced. 
For example, many authors have reported that the tran- 
sition temperatures of Al, In, Sn, and other soft metals 
can be increased over that of the bulk in the case of 
grains, films, or layered structures. [l8j 

Similar issues have been addressed theoretically in 
Hubbard models. For attractive models, spatial varia- 
tion in U can increase T c [l|,[2] for checkerboard and stripe 
patterns, especially when the modulation wavelength is 
close to the coherence length. [lj It is not surprising that 




(a) ID modulation 



(b)2D modulation 



FIG. 1: In this Letter we study the above inhomogeneity 
patterns, where certain bonds are made stronger (J s ) and 
others weaker (J w ) such that the total coupling is preserved. 



attractive-?/ Hubbard models benefit from inhomogene- 
ity. According to BCS theory, the pairing energy scale 
has a strong non-linear dependence on the attraction, 
A ~ T c ~ e~ x l vU . Since (e" 1 /^) > e^'"^, the lo- 
cal pairing amplifies favorable spatial variations in U . 
Even in the repulsive case, it has been shown that the 
superconducting gap of coupled 2-leg ladder systems is 
maximized for intermediate coupling between ladders. Q 
In this Letter, we predict ways to maximize the bulk 
transition temperature with respect to the zero temper- 
ature, long wavelength properties of the system. As a 
model for superconductors with low superfluid density, 
we study numerically 2D XY models with inhomoge- 
neous couplings. We are interested in patterns of the 
coupling constants that increase the transition tempera- 
ture over the homogeneous case. In order to make fair 
comparisons, we require that the zero temperature, long 
wavelength properties of the system remain unchanged. 
That is, we will not increase the low temperature energy 
density of the system, or the low temperature helicity 
modulus. For comparison, we also study Ising models 
with inhomogeneous couplings; the results support our 
findings for XY models. We find that although most 
patterns of inhomogeneity reduce the transition temper- 
ature T c , there are indeed certain "framework" patterns 
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of inhomogeneity that increase T c by up to a theoretical 
maximum of 76%. 




(a)Dimensionless helicity modulus u = T/T as a function 
of system size. The curves are solutions of the scaling 
equations, Eq. J3}, chosen to fit the data (crosses). The 
dashed line is T/T = — . 
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(b)Helicity modulus as a function of temperature. The 
black curve is T(L = oo,T) obtained by FSS. The dashed 
line is T = —T. 



FIG. 2: Finite-size scaling. The crosses are Monte Carlo re- 
sults for the helicity modulus T(L, T) of a 4x4-modulated XY 
model with J w — and J s — 4J avg , for T — 0.9, 1.0, 1.5 in 
units of J avg and L = 4,8, ...1024. 



The XY model has the following classical Hamiltonian: 

(1) 



n XY [e 



where i,j are site labels, Jij are nearest-neighbour cou- 
plings, and 9 X are real-valued phase (angle) variables. We 
choose the couplings Jy so as to preserve Ylij Jij- This 
leaves the zero-temperature properties of the system un- 
changed, such as the helicity modulus T(T = 0) and the 
energy U(T = 0). We consider here only two-dimensional 
models. 

We are interested in how to optimize the transition 
temperature, and to study this we focus on the behavior 
of the helicity modulus T(T) of the XY model, which is 
directly proportional to the superfluid stiffness n s of a 
phase-dominated superconductor, (lij by We study this 
quantity via Monte Carlo simulations on square lattices, 



where T(T) can be calculated using 



2V \ ^ 13 

\ (ij) 



1 2> 



^2 J ij sin (0i - 6 i) 

(ij) 



(2) 



We use the Wolff cluster algorithm [4(, which is the fastest 
serial algorithm for our purposes. The 9 variables are 
stored and manipulated as two-vectors to avoid trigono- 
metric function calls. 

In order to obtain reliable estimates of T c , we have 
performed finite-size scaling (FSS) on T(L,T) in the fol- 
lowing manner. The KT transition can be described by a 
two-parameter scaling flow 0, 0, 0, H| for the dimension- 
less helicity modulus u — T/T and the 'vortex fugacity' 

y, 



da, 
HI 



-4:TT 3 u 2 y 2 , 



^ = (2-™)y, (3) 



where I = hiL is the length scale. This pair of differential 
equations can be solved numerically, given initial values 
u(lo) = ito and y(Zo) = Vo (where l is some reference 
length scale). For each temperature T, we choose uq and 
yo so as to obtain a good fit of u(l) to the Monte Carlo 
data for the available system sizes, 4 < L < 1024 (see 
Fig. 2(a)). We then integrate the differential equations 
all the way to I — oo. This gives u(oo) and hence the 
helicity modulus in the infinite-size limit T(T,L = oo), 
shown in Fig. [2(b)] [H 




FIG. 3: Y(T, L = oo) for 2D XY models. The red curve is 
for the homogeneous 2D XY model. The green curve is for a 



4x4 modulation with J s tron 



3.4 and J w 



0.2. The blue 



curve is for a 4x4 modulation with J s tron g = 4 and J wca k = 0. 
The dashed line is T = -T; the arrow indicates a theoretical 
upper bound on T c , T^ax = f Javg- 

Our most important result is that by redistributing the 
bond strengths of an XY model in certain inhomogeneous 
patterns, it is possible to increase T c . As a concrete ex- 
ample of how this comes about, we show how the shape 
of the helicity modulus curve vs. temperature is changed 
by introducing inhomogeneity. In Fig. [31 we show our 
simulations of T(T), extrapolated to infinite system size 
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L — > oo, for a 2D inhomogeneity of the type shown in 
Fig. 1(b) using A — 4a where a is the underlying lattice 



like those in Fig. 1(a) Since this type of inhomogeneity 



constant. In the blue and green curves of Fig. [31 the cou- 
pling constant Jy has been made stronger on the darker 



lines in Fig. 1(b) (J s trong), and weaker on the lighter 
lines (J W eak)- We compare these to the uniform case (the 
red curve) with J set equal to the spatial average of J wca k 

and Jstrong, J — Javg = (Jstrong + (A — 1) Jwcak) /A. ThuS 

for all curves shown in Fig. [3l the zero temperature he- 
licity modulus Y(T = 0) and the zero temperature free 
energy are the same. 

In the uniform case, it is known that T c = 0.8929J 
[1, El , and that the low temperature slope of the helic- 
ity modulus T'(0) = 1/4 (Tfl E3- Thc 

green curve m 

Fig. [3] shows the helicity modulus for Jstrong = 3.4 and 
Jweak = 0.2. In this case, the transition temperature 
is enhanced by 8% above the homogeneous case. For 
the case of extreme inhomogeneity with Jstrong = 4 and 
Jweak = (the blue curve), the transition temperature 
is 25% higher than in thc uniform case. The shape of 
the green curve demonstrates the separation of energy 
scales that happens with inhomogeneity. Notice that at 
the very lowest temperatures, the green curve is dom- 
inated by the long-wavelength average of the coupling 
constants, and the low temperature linear slope of the 
helicity modulus is identical to that of the homogeneous 
case. As temperature is raised, the slope increases in 
magnitude, as the weak plaquettes become disordered. 
Then, at a higher temperature, the slope approaches that 
of the blue curve, indicating that at higher temperatures 
the helicity modulus is dominated by Jstrong- It is this 
shallower high temperature slope which causes the helic- 
ity modulus to overshoot the homogenous T c , and leads 
to an inhomogeneity-induced enhancement of the transi- 
tion temperature. 

For 2D patterns like those in Fig. |l(b)[ the enhance- 
ment of the transition temperature increases with A, as 
shown by the green curve of Fig. |4(b)| However, the 
enhancement is constrained by the zero temperature he- 
licity modulus. Even in the presence of inhomogeneity, in 
2D the system remains in the KT universality class, and 
the helicity modulus has a universal jump at the tran- 
sition such that T(T C ) = 0.6365T C . Since thermal fluc- 
tuations introduce disorder, Y(T C ) < T(T = 0), so that 
T c /T(0) < 1.57, or equivalents, T c /T c0 < 1.76. This the- 
oretical upper bound on T c is illustrated in Fig. [31 That 
is, although the zero temperature properties of the sys- 
tem may be used as a predictor of the transition temper- 
ature, T c = AT(T = 0), the constant A is non-universal. 
In fact, A = T c /T(0) may be useful for characterizing the 
degree of inhomogeneity: it increases from 0.89 for the 
uniform XY model up to a theoretical maximum of 1.57. 
For a 2D system, a large measured value of this ratio may 
indicate substantial inhomogeneity. (An increase in this 
ratio may also indicate higher dimensionality. [10j ) 

We have also considered one-dimensional modulations, 



drives the system towards more one-dimensional physics 
where a phase transition is forbidden by the Mermin- 
Wagner theorem, the transition temperature decreases, 
as shown in Fig. 4(a) Hence the enhancement of T c 
is not additive — thc effect of a 2D modulation is not 
double that of a ID modulation. 

Fig. [4] also shows the effect of inhomogeneity in the 
Ising model, for comparison. Inhomogcncous Ising mod- 
els may be described by the Hamiltonian 
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sing 



[a] = - Ji^OiO 

(ij) 



J > 



= ±1. 



(4) 



As with the XY model, we restrict ourselves to two di- 
mensions. For the purpose of studying different length 
scales of inhomogeneity, we focus on an extreme type 
of inhomogeneity with J wea k = 0, and Jstrong = AJ avg - 
Such patterns correspond to 'decorated' lattices. By inte- 
grating out all doubly-coordinated spins (that is, by ap- 
plying the so-called decorated-iteration transformation, 
J c ff = tanh _1 (tanh Ji tanh J2)), one can reduce a deco- 
rated lattice to a primitive lattice and thus obtain an 
exact expression for its T c (Eq. {5}). [21] (Unfortu- 
nately, the decoration-iteration transformation for XY 
models involves an infinite set of Fourier components of 
the potential [7( and it does not lead to exact results for 

In Fig. [4) we show the effect of extreme inhomogeneity 
(i.e., with Jwcak = 0) on the transition temperatures in 
Ising and XY models. We use the maximum value of 
Jstrong = AJavg, because for a given wavelength A, this 
gives the largest enhancement of T c while conserving the 
average coupling J avg - While we are interested primarily 
in superconductors with small superfluid density, which 
can be captured with an XY model, we also show results 
for the Ising model, for which results can be obtained 



analytically as described above. Fig. 4(a) shows the effect 



of a purely ID modulation, as a function of distance A 
between strong bonds Jstrong, chosen so as to preserve 
the zero temperature, long wavelength properties of the 
system. The pattern of coupling constants is shown in 
Fig. 1(a) In the Ising case, the transition temperature 
is unchanged by this procedure. In the XY case, the 
transition temperature decreases monotonically with A. 

The effect of a 2D modulation is shown in Fig. |4(b)[ 
Again, parameters are chosen so as to preserve the zero 
temperature properties of the system. Fig. l(b)| shows 
the pattern of coupling constants. Here, the transition 
temperature in the Ising case increases as 



2A J a 



sinh 1 [csch(i sinh 1 1)] 



(5) 



One of the occurrences of lambda in this equation is due 
to taking A bonds in parallel, to form "bundles" , and the 
other occurrence is due to taking A bundles in series. 
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For XY models, the transition temperature also increases 
monotonically with modulation length A. In this case 
there is an upper bound, set by the zero temperature 
properties of the system, as shown in Fig. [3] That is, the 
maximum enhancement of T c possible with this type of 
inhomogeneity in an XY model is 76%. 

In conclusion, we have shown that certain types of in- 
homogeneity can increase the transition temperature of 
Ising and XY models. Specifically, two-dimensional mod- 
ulations of the coupling constants that preserve the spa- 
tial average coupling increase the transition temperature 
over that of the uniform case. One-dimensional modu- 
lations depress the transition in XY models, and leave 
the transition temperature unchanged in Ising models. 
Our results for 2D XY models may indicate that certain 
types of inhomogeneity can result in an enhancement of 
superconductivity in systems with low supcrfluid density. 
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FIG. 4: Critical temperatures T c of the lattices depicted in 
1(a) and 1(b) as a function of the wavelength of the 
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inhomogeneity. 
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The helicity modulus measures the change in the free 
energy caused by a small change in the phase angle [!|. 
This treatment neglects the tiny corrections due to non- 
zero winding numbers 9] , but it is adequate for the level 
of accuracy of the present work. As a check, we have esti- 
mated T c using other methods, such as finite-size scaling 
for the susceptibility \{L, T) and for the Wolff cluster 
size N c i uat (L,T), based on the predictions of Kosterlitz- 
Thouless theory that x{L,T c ) ~ L 7/4 (lni) 1/s and that 

x(oo,T - T c ) ~ t- 1/16 e ci ~ 1/2 . The results agree with 
those obtained from finite-size scaling for T. 
For more complicated patterns not amenable to 
decoration-iteration, T c can be calculated using either the 
Pfafhan method 15] (e.g., for periodic inhomogeneity), or 
the recently developed bond-propagation algorithm [l6l] 
(which is the most efficient in the case of bond dilution) . 



